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Abstract 

The scaling form of the free-energy near a critical point allows for the definition 
of various thermodynamical amplitudes and the determination of their dependence 
on the microscopic non-universal scales. Universal quantities can be obtained by 
considering special combinations of the amplitudes. Together with the critical ex- 
ponents they characterize the universality classes and may be useful quantities for 
their experimental identification. We compute the universal amplitude ratios for 
the Tricritical Ising Model in two dimensions by using several theoretical methods 
from Perturbed Conformal Field Theory and Scattering Integrable Quantum Field 
Theory. The theoretical approaches are further supported and integrated by results 
coming from a numerical determination of the energy eigenvalues and eigenvectors 
of the off-critical systems in an infinite cylinder. 



PACS numbers: 64.60.Fr, 05.50+q,75.10Hk 



1 Introduction 



One of the most powerful and fascinating concepts in the investigation of critical phe- 
nomena - which has successfully passed the scrutiny of both experimental and theoretical 
tests during the last decades - goes under the name of universality According to this 
principle two statistical models, which share the same symmetry of the order parameters 
and the dimensionality of the space of their definition, show an identical critical behavior 
although they may greatly differ in their microscopic realizations: near the phase transi- 
tion, when the correlation length is much larger than any other microscopic scale, they 
appear as two representatives of the same universality class. The first characteristic of 
a given universality class consists of a set of critical exponents. Their values are gener- 
ally given in terms of algebraic expressions of the anomalous dimensions of the relevant 
operators present at the critical point. In two dimensions, by the powerful methods of 
Conformal Field Theory (CFT) P, |], ^, |^ one can ascertain both the values of the critical 
exponents, the Operator Product Expansions (OPE) and the multi-point correlators of 
the relative fields. 

However, a complete analysis of the class of universality should also include the de- 
scription of the structure of the Renormalization Group trajectories near the critical point. 
The most ambitious goal would be the determination of both the scaling function which 
fixes the equation of state and the off-critical correlators of the various order parameters^. 
Although the exact determination of the equation of state of a given universality class may 
often be a difficult task, the scaling property alone of the free-energy is nevertheless suf- 
ficient to extract numerous predictions on universal combinations of critical amplitudes. 
As it will become clear in section 3, these universal combinations are pure numbers which 
can be extremely useful for the identification of the universality classes. In fact, the am- 
plitude ratios are numbers which typically present significant variations between different 
classes of universality, whereas the critical exponents usually assume small values which 
only vary by a small percent by changing the universality classes. Hence the universal 
ratios may be ideal experimental marks of the critical scaling regime 0. 

In recent years, due to the theoretical progress achieved in the study of two-dimensional 
models, some universal ratios and other universal quantities have been computed for a 
large variety of bidimensional systems, such as the self-avoiding walks P], the Ising model 
[9-14] or the q-state Potts model ||15| , [1^, to name few. In this paper we will concen- 
trate our attention on the determination of the amplest possible set of universal amplitude 
ratios relative to the class of universality of the two-dimensional Tricritical Ising Model 



^Near the critical point all quantities relative to different models of the same class of universality 
become identical provided an opportune rescaling of the order parameters, the external fields and the 
correlation length of the models is made. 
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(TIM). Preliminary results relative to some off-critical phases of this class of universality 



have been presented in our previous publication and the aim of this paper is twofold. 



First of all, to complete the list of universal amplitude ratios of the TIM presented in |18 



and to perform an exhaustive analysis of all its possible phases. Secondly, to illustrate 
in full detail the theoretical methods which have been employed in such determination: 
in view of their successful applications, these techniques may be useful to analyze and to 
obtain similar results for other statistical models, in such a way to bridge a closer contact 
between theoretical and experimental results in two-dimensional physics. 

Introducing the Tricritical Ising Model, from a field theoretical point of view this 
model may be regarded as a Landau-Ginzburg (LG) $^-theory near its tricritical point 
1^. The LG terminology allows a qualitative understanding of the phase-structure of 



the model effortlessly. However the LG approach is often too elementary for the under- 
standing of some remarkable symmetries present in the 2-D TIM. As a matter of fact, 
the bidimensional TIM is an unique example of critical phenomena: it is still sufficiently 
simple to be solved but at the same time it presents an extremely rich and fascinating 
structure of excitations which can attract the curiosity of a theorist. Depending on the 
direction in the phase space in which the system is moved away from criticality, one can 



observe, for instance, a behavior ruled by the exceptional root system E-j [g^j HI' ^ 



or by supersymmetry [^, ^ ^ (in its exact or broken phase realization) or by an 
asymmetrical pair of kinks |28|, |29| . In addition, the description of its low-temperature 
phase is easily obtained from the one of its high-temperature phase because of the self- 
duality of the model. From the experimental point of view, a number of physical systems 
exhibit a tricritical Ising behavior, among them fluid mixtures or metamagnetsQ. Hence, 
there is an obvious interest in computing the ampler set of data for this class of univer- 
sality and in testing the theoretical predictions versus their experimental determination. 
Our calculation of the universal ratios of the TIM will be performed by a combined use 
of results coming from Perturbed Conformal Field Theory, from the integrable structure 
of some of the deformations of the critical point action and also from some numerical 
approaches. 

The paper is organized as follows. In Section 2 we briefly describe the universality class 
of the TIM and the symmetry properties of the theories resulting from the deformations 
of the critical point action made of each of the four relevant fields. Section 3 is devoted 
to the discussion of the scaling behavior of the singular part of the free-energy and the 
definition of the universal ratios obtained by considering some particular combinations 
of the thermodynamical amphtudes. In Section 4 we discuss the Quantum Field Theory 
approach to the computation of the universal ratios. The numerical method based on the 

^The interest reader may find ample review of this topics in rcf. PQ|. 
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diagonalization of the off-critical hamiltonian obtained by truncating the Hilbert space 
of the conformal space is discussed in Section 5. A thorough analysis of each relevant 
perturbation of the TIM is performed in Section 6. Finally our conclusions are presented 
in Section 7. Several appendices, relative to some technical aspects of our calculations, 
are included at the end of the paper. 



2 The Class of Universality of the Tricritical Ising 
Model 

In this section we will briefly outline the main properties of the universality class of the 
two-dimensional Tricritical Ising model, whereby the detailed discussion of its physical 
properties relative to each of its perturbations is contained in Section 6. 

A possible lattice realization of the Tricritical Ising model is provided by the so-called 



Blume-Capel model This involves two statistical variables at each lattice site, Sk - 
the spin variable - which assumes values ±1 and tk - the vacancy variable - with values 
or 1, which specifies therefore if the site is empty or occupied. It is characterized by the 
most general Hamiltonian with nearest neighbor pair interaction 

N N N N N 

n = -Jj2 SiSjUt, +AY,ti-HY^ SiU - Hs ^{sititj + Sjtjti) -KJ2 Utj ■ (2.1) 

The parameter H represents an external magnetic field, an additional subleading 
magnetic source, J the coupling between two nearest occupied sites, A the chemical 
potential coupled to the vacancies and K an additional subleading energy term between 
them. Another possible two-dimensional lattice realization of this class of universality is 
provided by the so-called dilute models, discussed in |^ . 



By adopting a field theoretical point of view, a convenient way to analyze the univer- 
sality behavior of the TIM consists in considering a Landau-Ginzburg formulation based 
on a scalar field $(x) [0]. The main advantage of this approach is an account of the Z2 
symmetry properties of each order parameter which provides an easy way of understand- 
ing the phase structure of the model, at least qualitatively. In this formulation, the class 
of universality of the TIM is associated to the Euclidean action 

■1 



A= I d^x 



(2.2) 



with the tricritical point identified by the bare conditions 9i = 92 = 93 = 94, = 0. In a 
close comparison with the Blume-Capel lattice formulation of the model, the statistical 
interpretation of the coupling constants is as follows: 91 plays the role of an external 
magnetic field h, 92 measures the displacement of the temperature from its critical value 



3 



(T — Tc), may be regarded as a staggered magnetic field h' and finally may be thought 
chemical potential for the vacancy density. 
A dimensional analysis shows that the upper critical dimension of the above LG model 
is D = 3, where the tricritical exponents are expected to take their classical values, 
excluding logarithmic corrections. The mean field solution of the model easily shows that 
the LG action (|2.2| ) has a tricritical point, i.e. a critical point where a line of a second 
order phase transition meets a line of a first order phase transition. Consider in fact the 
case where all the Z2 odd couplings are equally set to zero. The potential in this subspace 
is given by 

= ^2$^ + ^4$^ + '^>^ • (2.3) 
The line of a second order phase transition is identified by the condition 

92 = Q , (74 > , (2.4) 

whereas the line of a first order phase transition (where three degenerate vacua coexist) 
is obtained by the condition 

^2 > , (74 = -2^2 . (2.5) 

Hence the point (71 = (72 = fl'a = ^74 = is indeed a tricritical point. 

In two dimensions - the case which will mostly concern us - there are strong fluctua- 
tions of the order parameters and therefore the exponents and the amplitudes extracted 
by its mean field solution cannot be trusted. However, in this case one can take advantage 
of the powerful methods of Conformal Field Theory to obtain an exact solution of this 
model at criticality. In fact, the bidimensional TIM is described by the second represen- 
tative of the unitary series of minimal models of CFT p|, ^ |]: its central charge is equal 
to c = ^ and the Kac-table of the exact conformal weights of the scaling fields 

§5 ' l<k<A ^^-^^ 

is given in Table 1. There are six primary scalar fields 4>^-^, which close an algebra under 
the Operator Product Expansion 

<Pi{zi,Zi) <Pj{z2,Z2) ~ E C^ifc \Z1-Z2 |-2(A>+A,-A,) 0^(^,2,^2) • (2.7) 

k 

The skeleton form of this OPE algebra and the relative structure constants of the Fusion 
Rules of the TIM are in Table 2. The six primary fields can be identified with the normal 
ordered composite LG fields (see Table 3). With respect to their properties under 
the Z2 spin-reversal transformation Q : $ ^ — $ we have: 
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1. two odd fields: the leadiiiEf maErnetization operator a = (j)3_ j. = $ and the sublead- 

80 ' 80 

ing magnetization operator cr' = (f)j_ j_ =: : 

16 ' 16 

2. four even fields: the identity operator 1 = 0o,O; the leading energy density e = 
(hi 1 =: $^ :, the subleading energy density t = 6 & & =: $^ :, which in metam- 
agnets assumes the meaning of the density of the annealed vacancies, and the field 
e" =033. The OPE of the even fields form a subalgebra of the Fusion Rules. 

2 ' 2 

In the TIM there is another Z2 transformation - the Kramers- Wannier duality D - under 
which the fields transform as follows: 

• the order magnetization operators are mapped onto their corresponding disorder 
operators 

u = D-^aD = A , /i' = D-^a'D = d)^ jl ■ (2.8) 

80 ' 80 16 ' 16 

• the even fields are mapped onto themselves, 

D-^eD = -e , D-HD = t , D^V'D = -e", (2.9) 
i. e. e and e" are odd under this transformation whereas t is even. 

Interestingly enough, at criticality the TIM also provides a concrete realization of a 
supersymmetric theory since it is the first representative of the superconformal minimal 



models |^, ^: the even fields can be grouped into a superfield of the Neveu-Schwartz 
sector 

^^{z, z, e, e) = e{z, z) + e i){z, z) + e ip{z, z) + ee t{z, z) , (2.10) 

(where 9 and 9 are Grassman variables) while the magnetic fields give rise to two irre- 
ducible representations in the Ramond sectoiQ. The critical superconformal LG action is 
given by 

(2.11) 



A = j d^x(f9 
with the covariant derivatives defined as 



d d d — d 

At the critical point, the TIM can be also realized in terms of a coset construction 
of a Wess-Zumino-Witten model on the group G/H given by (-^7)1 ® {Ej)i/ {Ej)2 (for 
a general discussion on the coset model see, for instance [ffl). For the central charge 



^The disorder fields /i, /i' and the fcrmionic fields V'j "0 enter the partition function of the model on 
a torus with twisted boundary conditions, see for instance |^ where the relative Fusion Rules are also 
presented. 
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c = cg — cjy we have in fact c = 2x yjg — -^j = jq- Concerning the irreducible 
representations, for the WZW based on (-£^7)1 we have {I, Hg}i with conformal dimensions 
{0, |}, whereas for the WZW model (£'7)2 we have {I, Hi, lis, rig, n2}2 with conformal 
dimensions {0,^,|,|5,f|}- Therefore, the conformal fields of the TIM emerge from the 
decomposition 

X (i)i = (/)„M®(i)2+ f4^) ®(ni)2+f4) ®(n5)2 

X (He)! = (^) ® (n7)2 + ® (n6)2 (2.13) 

\lb/TiM \o(jJtim 

(He)! X (He)! = (|) x (1)2 

\Z/ TIM 

As will be discussed later, the above symmetries present at the critical point of the 
TIM are also useful for the investigation of some off-critical phases of the model. The 
four fields a, e, a' and t of increasing anomalous dimensions are, from a Renormalization 
Group point of view, all relevant operators (i.e. their conformal weight satisfies A < 1) 
and therefore they can be used to move the TIM away from criticality. To simplify the 
formulae below, it is convenient to adopt the compact notation [i = 1, 2, 3, 4) to denote 
collectively all these fields, so that (pi = a , = y^s = o'' and = t. In the vicinity of 
the critical point the partition function of the model can be expressed by the path integral 

^[^1,^2,^3,^4] = / P0e-[-^^^^+Sti9»/'^»(-)'^'-] = e-^~(^i'^^'^^'^^) . (2.14) 

We use the notation Zi[gi] = Z[gi, 0,0,0], Z2[g2] = Z[0, g2,0,0], etc. for the partition 
functions corresponding to the individual deformations of the conformal action. An imme- 
diate result for the off-critical phases can be drawn from the symmetry properties of the 
fields In fact, since the fields cpi and (p^ are odd under the spin-reversal Z2 symmetry, 
a change of the sign of the corresponding coupling constants gives rise to identical physical 
situations, i.e. Zi[gi] = Zi[—gi] and 2'3[5f3] = Z^l—g^]. The operator 02 is odd under the 
Z2 duality transformation and therefore the physical situations which originate from a 
change of sign of this coupling constant will be related by a duality transformation. Fi- 
nally the operator 04 is even under both Z2 symmetries and therefore the changing of the 
sign of the corresponding coupling constant will produce two distinct physical situations. 

At this stage it is also useful to anticipate the nature of Quantum Field Theories 
which originate from each individual deformation, postponing their detailed discussion in 
Section 6. The QFT associated to the deformation of the field (fi alone is not integrable: 
numerical indications of this fact were discussed in The QFT which originates from 
the deformation of ip2 is instead integrable and the pattern of its bound states and the 



scattering amplitudes are related to the hidden Ej algebraic structure of the model 



The deformation of the critical action by means of the field 993 produces an integrable 
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field theory made of kinks, whicli interpolate between two asymmetric vacua 
Finally, the <y94 deformation made with a positive value of the relative coupling constant 
corresponds to an integrable massless RG flow between the TIM and the standard Ising 
model ^ whereas the 994 deformation with a negative value of gives rise to an 
integrable massive QFT with kink excitations interpolating between three degenerate 



vacua pB]. Both these last QFT's give explicit realizations of a supersymmetric system 
in its broken and unbroken phase respectively. A useful summary of the theories resulting 
from each deformation can be found in Table 4. In conclusion, excluding the magnetic ipi 
deformation, all the others gives rise to integrable QFT. This fact will be quite important 
for our future considerations. 



3 Scaling Form of the Free— Energy and Universal Ra- 
tios 

The scaling property of the relevant fields which span the scaling region of the TIM near 
its critical point allows the derivation of a large set of universal quantities which are 
of experimental interest. For their derivation we will consider the general case of the 
Tricritical Ising Model defined in a D dimensional space, even though our final attention 
will be focalized on the two-dimensional system. At this stage of the discussion we do 
not take into account the eventual logarithmic corrections explicitly present in the 3- 
dimensional version of this modelQ and the possible ultraviolet renormalization effects. 
These, however, will be considered within the context of Sections 4 and 6. 

The scaling property of the order parameters is encoded into the asymptotic form of 
their two-point function^ 

(^,(x)^,(0)) ^ , |a;|^0 , (3.1) 

which therefore identifies the parameters Aj as the conformal dimensions of the fields. 
The standard conformal normalization of the fields is obtained by the choice Ai = 1. 
From the power law behavior of (|3.1| ) it follows that the coupling constants Qi behave as 

gi ~ A^-2^« , (3.2) 

where A is a mass scale. Therefore, moving the system away from criticality by means of 
one of the relevant field ipi, there will be generally a finite correlation length ^ which in 

^On this issue the interested reader may consult for instance js^ and references therein. 
^This notation is not standard for a generic D dimensional system but it has the advantage of an easy 
comparison with formulae which are valid in 2-D systems. 
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the thermodynamical limit scales as 

e~a(ir,(7,)"^ , (3.3) 

where a ~ may be regarded as a microscopic length scale. The terms Ki are metric, 
non-universal factors which depend on the unit chosen for measuring the external sources 
Qi, alias on the particular realization selected for representing the universality class. In 
the presence of several deformations of the conformal action, the most general expression 
for the scaling form of the correlation length may be written as 



{ = «...(/f,9.)---I£,^(^j , (3.4) 

where 

are the so-called crossover exponents whereas Ci are universal homogeneous scaling func- 
tions of the ratios There are of course several (but equivalent) ways of writing 
these scaling forms, depending on which coupling constant is selected as a prefactor. In 
the limit where gi Q {I ^ i) but gi ^ 0, equation (|3.4| ) becomes 



i^ = a^g, , "'^"^ . (3.6) 

Consider now the free-energy f[gi, . . . , g^. This is a dimensionless quantity defined 

by 

Assuming the validity of the hyperscaling hypothesis, in the thermodynamical limit its 
singular part (per unit of volume) will be proportional to the D power of the correlation 
length. Let us denote the singular part of the free-energy for unit volume by /[^^i, . . . , (74]. 
Depending on which scaling form is adopted for the correlation length, we have corre- 
spondingly several (but equivalent) ways of parameterizing this quantity 

/[^i, • • • ,^4] = /a^i, • • • , <74] ^ {K,g,)T^^ ( (/^y,, ) • (3-8) 

The functions Ti are universal homogeneous scaling functions of the ratios ^'^L-. . As 
it will soon become clear, there is an obvious advantage in dealing with different but 
equivalent expressions for the free-energy: in fact, since we will be mostly concerned with 
the physical situations which originate from pure deformations [i.e. those obtained by 
keeping only one coupling constant finally different from zero), the choice of which one 
has to be selected naturally follows from the particular deformation which is considered. 
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Let us discuss now the definition of the thermodynamical quantities related to the 
various derivatives of the free-energy. We will adopt the notation to denote expecta- 
tion values computed in the off-critical theory obtained by keeping (at the end) only the 
coupling constant gi different from zero. The first quantities to consider are the vacuum 
expectation values (VEV) of the fields ipj which can be parameterized as 

/ \ 



^B.^gr'^^ , (3.9) 

91=0 



with 



B.^^KjKr''^^ . (3.10) 
The above relations can be equivalently expressed as 



gi = Dij{{cpj),) , (3.11) 

with ^ 

Dij ~ c^riA- ■ (3-12) 



The generalized susceptibihties of the model are defined by 



dgk dgkdgj 



(3.13) 



91=0 



They are obviously symmetrical in the two lower indices. By extracting their dependence 
on the coupling constant gi, they can be expressed as 

^k-^k9^ ""'"^ , (3-14) 

with 



r),^K,K,K, "-'-^ . (3.15) 

Some of the above quantities have, of course, a very familiar meaning. For instance {(fii)i 
is nothing but the mean value of the magnetization in the off-critical theory defined by 
the i-th deformation while F^.! is the associated magnetic susceptibility. Similarly, {(f2)i 
is the mean value of the energy along the i-th deformation of the critical theory and 1^22 
the specific heat. 

As easily seen from the above formulae, the various quantities obtained by taking 
the derivatives of the free-energy contain metric factors (the quantities Ki) which make 
their values not universal. However, it is always possible to consider special combinations 
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thereof in such a way that all metric factors cancel out. Here we propose the consideration 
of the following universal ratios 

{Rc)], = ; (3.16) 
^| = (ry'^''e° ; (3.18) 

{Ra)) = rj, A. ■ (3.19) 

= rf (|j • (^-'^^ 

These quantities are pure numbers which therefore characterize the universality class of 
the model. Their definitions closely follow and generalize the ones relative to the fa- 
miliar Ising model (see, for instance P, |12|). Other universal ratios may be defined as 
well and in fact some of them will be considered in Section 6 devoted to the analysis of 
each deformation of the critical point action. Since we will individually compute all the 
important quantities involved {Bij, F*^, etc.), there is really no problem in considering 
other universal combinations, if one wishes to do so. It is worth emphasizing that, from 
an experimental point of view, it should be simpler to measure universal amplitude ratios 
rather than critical exponents: in fact to determine the former quantities one needs to 
perform several measurements at a single, fixed value of the coupling which drives the 
system away from criticality whereas to determine the latter, one needs to make mea- 
surements over several decades along the axes of the off-critical couplings. Moreover, 
although not all of them are independent, the universal ratios are a larger set of numbers 
than the critical exponents and therefore permit a more precise determination of the class 
of universality. 



4 Quantum Field Theory Approach 



Essential quantities of the universal amplitude ratios ( p.l6| ) - (|3.20| ) are the correlation 
length prefactor the VEV amplitudes Bji and the generalized susceptibilities Vjj^. 
With the aim of determining these quantities, in this section we will discuss some useful 
results relative to the two-dimensional Quantum Field Theories (QFT) associated to the 
Renormalization Group flows originating by relevant deformations of the conformal action. 
These QFT are, after all, particular representatives of the universality class of the model 
and from now on we will only focalise on them for the study of the off-critical dynamics. 
The advantages of adopting this approach will soon become evident. 
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In the following we assume the fields to be normalized according to the conformal 
normalization, i.e. 



'■5 

|4A 



lim Jxl"^' ((^i(x)(^,(0)) = 1 , (4.1) 
and we denote by Aif the QFT associated to the action 

A = AcFT±9i J Viix)d^x , gi>Q . (4.2) 

The coupling constant Qi is a dimensional quantity which can be related to the lowest 
mass-gap mj = ^^"^ of the theory according to the formula 

^7, = C,m^2^' , (4.3) 

or, equivalently 



mi = C,gr^ , (4.4) 



with Cj = Cj ^ . When the QFT associated to the action (f4.2|) is integrable, the 
pure number Cj can be exactly determined by means of the Thermodynamical Bethe 



Ansatz ||3^, |3^. When the theory is not integrable (this is the case for the magnetic 
deformation of the TIM), the constant Ci can be nevertheless determined by a numerical 
method, based on the so-called Truncated Conformal Space Approach p8|, which will be 



discussed in Section 5. In conclusion, for all individual deformations of the TIM we are 
able to completely set the relationship which links the coupling constant to the mass-gap 
of the theory and therefore we are able to switch freely between these two variables. 

Another set of quantities which can be fixed by QFT are the matrix elements of the 
order parameters, the simplest ones being the vacuum expectation values. In this case we 
have 

{<^j)i = Bji mf^' , (4.5) 



I.e. 



{^jh = B,igr^^ , i?,, = S,,cf^' . (4.6) 



When the theory is integrable, the constant Bji can be fixed exactly, thanks to the results 
of a remarkable series of papers 0. When it is not integrable, the constant Bji can 



be nevertheless estimated by means of a numerical approach, first proposed in [^, which 



will be reviewed in Section 5. Hence, also in this case, we are able to determine completely 
these quantities. Let us present the exact expressions of the VEV for the three integrable 



deformations of the TIM, which are obtained by specializing the formulae of ref. In 
the expressions below, the fields are labelled by their position {I, k) in the Kac table of 
the model (see eq. ( p.6| ) and Table 1) and the lowest mass-gap of the different theories is 
simply denoted by m 
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• For the (/72 energy deformation we have 



(0. I ^i,k I 0,) = 



sin [f\5l- Ak l) 



sin^ 



5 mTrF 



2A, 



_2!v^r(i)r(|) 

where, for \ Rer] \< 4, Qi,2(^) is given by the integral 

^ , , ( dt f sinhGf siiili(f(// - 1)) siiilifff// + 1)) 



Qi,2{5l-ik) , (4.7) 



sinh 15t sinh lOi sinh At 



X 



f„2 _ i\ N 

(cosh 18t + cosh 8t - cosh 16t + cosh 4t + 1) - ^ ^ ^ e~^* 



(4. 



and is defined by its analytic continuation outside that domain. The index s labels 
the various vacua and takes different values depending on the sign of the coupling 
constant for g2 < the spin symmetry a — > —a is spontaneously broken and 
there are two vacua identified by s — 1,3; ior g2 > there is a unique ground state, 
associated to s = 2. 



For the 993 sub-leading magnetic deformation we have 



(0. I I 0, 



sm(^ I r>l - ik 



sm 



1 in T 


T 


(5) 






(1) 


r 


(1) 



2A 



l,k 



Q2,i(5/-4A;) , (4.9) 



where, for | iier/ |< 5, Q2,i{v) is given by the integral 



^ , , ( f°° dt f smh3t smh(t(r7 - 1)) smh(t(r7 + 1)) 

Q2i(^)=exp</ — — —X 

'^'^ ^[Jo t \ sinh 9i sinh 5t sinh 8t 

(cosh 9t + cosh t - cosh lit + cosh 5t + l)- ^ ^^ e"^*j | , (4.10) 

and is defined by its analytic continuation outside that domain. For this deformation 
there are two vacua associated to the values s = 2, 4. 

• For the </?4 vacancy density deformation, in its massive phase we have 



(0. I ^i,k I 0,) 



sin {f\5l- Ak 1) 



sinf 



m0Fr(|) 



2A, 



Qi,3(5/-4A;) , (4.11) 



where, for | Re 7] |< 4, Qifi{ri) is given by the integral 



f rft /cosh2t sinh(t(?7 - 1)) sinh(t(r7 + 1)) {rf - 1) 
Qi,3(77) - exp 2 cosh t sinh 4t sinh5t 40~' )] 

(4.12) 

and is defined by its analytic continuation outside that domain. For the massive 
phase of this deformation we have three vacua labelled by s — 1, 2, 3. 
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As discussed in Section 5, in addition to the above vacuum expectation values, a 



generalization of the numerical approach of ref. |^T[] often leads to a reasonable estimation 
of the matrix elements of the order parameters between the vacuum states and some of 
the excited states, as for instance {0\ipj\Ak)i where Ak is a one-particle state of mass M^. 
These quantities will be useful for obtaining sensible approximation of the large-distance 
behavior of several correlators. 

Another useful piece of information on the off-critical dynamics can be obtained by 
exploiting the properties of the stress-energy tensor T^i^{x). In the presence of the per- 
turbing field (fi, the trace of the stress-energy tensor is different from zero and can be 
expressed as 

e{x) =2ng,i2-2A,)ip, . (4.13) 
The vacuum expectation value of 0(x) is given by 



{e) = w^mj = w,gr^' , (4.14) 

with Wi = WiCf. As before, when the theory is integrable the constant Wi can be deter- 
mined exactly, otherwise it can be computed numerically]^. The trace of the stress-energy 
tensor enters two useful sum rules which link conformal data of the ultraviolet fixed point 
to off-critical quantities. The first of them - called the c-theorem sum rule ||^ - re- 
lates the central charge c of the ultraviolet theory to the second moment of the two-point 
connected correlation function of G 

c = l-Jd\ ix^eix) e(o))e . (4.15) 

For all relevant deformations, it is easy to check that the above integral is always conver- 



gent. The second sum rule - called the A-theorem sum rule |43| - reads 



i.e. it relates the conformal dimension Aj of the field (pj to its VEV and to the integral of 
its connected off-critical correlator with G(x). This time the above integral is not always 
convergent (the detailed analysis of its convergence may be found in the original paper 
[^ ). Notice, however, that the A-theorem also involves the VEV of the field (pj and it is 
easy to see that the divergence/convergence of the integral is always accompanied by the 
divergence/convergence of the VEV in such a way that the sum-rule always mantains its 
validity. The proof of this statement is simple: by taking the derivative of the VEV (|4.6| ) 
with respect to gi we have in fact 



^Obviously this number can be also obtained in terms of the VEV of the field (pi and the expression 
of given in eq. (4.13). 
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On the other hand, the above quantity can also be computed by means of the fluctuation- 
dissipation theorem and is given by 

— (y.,), = - / d'x{^,ix)^,mi . (4.18) 



dgi 

By using eq. ( |4.13| ) and comparing the two expressions, the divergent /convergent nature of 
the integral is therefore directly linked to the divergent /convergent nature of the constant 
Bji entering the VEV of the field (pj. These considerations suggest that the quantity Aj 
on the left hand side of ( |4.16| ) is in any case obtained, also when the integral of the two- 
point function (B(x)<y9j(0))* diverges. In this case one needs to perform a simultaneous 
analytic continuation of both the integral and the corresponding VEV. 

Basic quantities in the universal ratios are the generalized susceptibilities F*^. By 
using equations. ( pl7|) and ( p.l4|) , the fluctuation-dissipation theorem provides the relation 
between f and the integral of the connected correlator 

t% = I £x{^^{x)MO)Yc ■ (4.19) 

The dependence on the coupling constant Qi of these quantities can be easily extracted. 
In fact, the connected correlator can be parameterized as 

(^.(^)V^.(O))^. = ^:^KJT^QU^r) . (4.20) 
(r =1 X I). Its dependence on m is obtained with a change of variable and by using the 
relation ( [4.4|) , we finally have f*;, = Vji^g^ with 

= Cf ^-^^^^-^ / dr-,^Q,,{r) . (4.21) 

Some of the above susceptibilities can be determined exactly, such as the components F-^, 
whose values are provided by the A-theorem sum rule 

= -T^B,, . (4.22) 

In all other cases, when an exact formula is not available, our strategy to evaluate the 
generalized susceptibilities will rely on two different representations of the correlators. 
These representations have the advantage to converge very fast in two distinct regions: 
the first representation is based on Conformal Perturbation Theory and allows a very 
efficient estimation of the correlation function in its short distance regime, while the second 
representation is based on the Form Factors and allows an efficient control of its large 
distance behavior. Due to the fast convergent nature of the two series in their respective 
domains, they are efficiently approximated by their lowest terms, which therefore can be 
evaluated with a relatively little analytical effort. These considerations obviously lead to 
the estimation of the integral ( |4.19| ) according to the following steps: 
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1. Express the integral in polar coordinates as 

r;., = 27ry^ rfrr(y,,(rK(0))^ , (4.23) 
and split the radial integral into two pieces as 

r+oo rR r+oo 

1= drr{^,{r)Mm = / drr{^,{r)MO))[+ dr r {^,{r)MO))l 

Jo Jo J R 

= hiR) + l2{R) . (4.24) 

2. Use the best available short-distance representation of the correlator to evaluate 
Ii{R) as well as the best available estimate of its large-distance representation to 
evaluate /2(-R). 

3. Optimize the choice of the parameter R in such a way to obtain the best evaluation 
of the whole integral. In practice, this means looking at that value of R for which 
a plateau is obtained for the sum of Ii{R) and I2{R)- 

In order to proceed in the above program it is useful to briefly recall the main features of 
the short-distance and long-distance expansions of the two-point correlation functions. 



4.1 Short— Distance Expansion 

A clear discussion on the perturbative ultraviolet renormalization of the fields and on 
the short-distance expansion of the two-point functions can be found in the references 



Here we will briefly review the main results useful for our purposes. 
First of all, there is a one-to-one correspondence between the fields at and away from 
criticality. However, renormalization effects induced by ultraviolet divergences can have 
the effect of expressing the off-critical fields in terms of a combination of the critical ones. 
Let us denote by $i(x) and $i(x) the conformal and the off-critical fields respectively. 
Consider the off-critical action obtained by a perturbation of a relevant (scalar) conformal 
field (A^ < 1) 

A = AcFT + 9 J d^x ^{x) . (4.25) 

In a conformal perturbative evaluation of the correlators which involves one of the field 
$i(0) we have 

(...$,(0)) = {..AmcFT + 9 f d^x{..MOM^))cFT + ... (4.26) 

Je<\x\<R 

where e and R are the ultraviolet and the infrared cutoffs. Let us analyze first the 
ultraviolet behavior of the above integral. This is controlled by the OPE of the two 
conformal fields 

4>(a;)<i>,(0) = E ^1^ I ^ |2(A,-A,-A,^) ^^(Q) ^ ^4 27) 
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and therefore the integral in ( 4.2(: ) is divergent if in the above expansion there are fields 
Afc(x) whose conformal dimensions satisfy the condition 



7fe = Afc - - A^^ + 1 < . (4.28) 

If this is the case, the off-critical renormalized field which has a finite correlator at the 
lowest order in g is defined by 

= (7vr 5: ^ 6^^'= A, + 0{g') . (4.29) 

k 7fc 

Hence, due to the ultraviolet divergences there may be a mixing of the initial conformal 
operators with a finite numbers of fields of lower conformal dimensions. 

The conformal perturbation series for the correlation functions also suffers from in- 
frared divergences. These divergences, however, cannot be absorbed into a redefinition 
of the local fields and, as a result, we have a non-analytic dependence on the coupling 
constant g. This non-analytic behavior is essentially due to the non-adiabatic change 
of the vacuum state in passing from the Conformal Field Theory of the fixed point to 
the generic massive theory of the off-critical system. To overcome this difficulty, one can 
adopt the strategy of considering the off-critical OPE 

^,{x)y^,{0) = J2Cf^{g;x)ApiO) , (4.30) 
p 

where the Ap{x) belongs to a complete set of local fields of the theory {i.e. the perturbed 
version of the conformal fields Ap{x)) whereas the structure constants Cfj{g; x) are analytic 
in g (as expected by their local nature). In this way all the non-analytic behavior of the 
correlation function 

(^,(x)^,(0)) = j:CUg;x){Ap{0)) (4.31) 

i 

is completely encoded inside the non-perturbative VEV's {Ap{0)) = Apg^^^^' . Concern- 
ing the structure constants Cji.{g;x), by dimensional reasons they admit the expansion 

oo 

CUg; x) = r2(^-^-^.) Y: Clf{gr'-'^-r , (4.32) 

n=0 

where r =| x | and they can be computed perturbativelyQ in g. Their first order contri- 
bution is given by |^ 

Cif = -fd^w {AP{ooMw)ip,il)ip,iO))cFT , (4.33) 

''The above formula has to be opportunely corrected when there is a resonance phenomenon among 
the fields of the conformal families. 
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where the prime indicates a suitable infrared (large distance) regularization of the inte- 
gral. It can be calculated by means of different approaches (as, for instance, a minimal 
subtraction scheme based on the OPE |^5| or an analytic prolongation in the parameters 



of the integrand of ( [4.33| ), i.e. the conformal weights). As shown in an efficient way to 
extract the finite part of the integral is provided by the Mellin transformation, discussed 
in Appendix A. This is the approach which we have mostly used in our calculations. 

In conclusion, by using the VEV of the fields and the first approximation of the 
structure constants of the OPE, we can obtain a reasonable approximation for the short- 
distance behavior of the connected two-point functions entering the definition of the 
susceptibilities in terms of the expression 

(^.(-)V^.(O))^ = T.t X'^Zi.) ^t\9r'-'''r + 

p n=0 

-{Vt)k{Vj)k ■ (4.34) 

The index k indicates which perturbation is considered. Since the short -distance rep- 
resentation ( [4. 311) is an expansion in the parameter a truncated form of the series 
( [4.34| ) is expected to be sufficiently accurate for r <^ ^. However, the convergence of the 



truncated series is often much better and results being sufficiently accurate also for r ~ 
as confirmed in several examples (see, for instance [|lT], [l^, 0). Hence, the above 

truncated form ( |4.34| ) can be confidently used for the evaluation of the integral Ii{R) 



in eq. ( 4.24 ). Finally, notice that the short-distance expansion of the correlators can be 
implemented independently on the integrable or non-integrable nature of the off-critical 
theory. 

4.2 Large Distance Expansion 

An efficient way to control the behavior of the correlators in the opposite regime (^jj ^ 1 
is provided by their spectral representation expansions. In this approach, one makes use 
of the knowledge of the off-critical mass spectrum of the theory to express the correlators 

{^i{x)^,m = Y.gn{r) , (4.35) 

n=0 

where 

9n{r) = I ^...^(O|y.,(O)|A„,(0O---^a„(^n))x (4.36) 

x(A,,(^i) . . .A„„ (e„) 1^.(0) |0) e-^EL.-'^^^he, _ 



The expression (4.36) has to be suitably modified in presence of kink excitations or massless particles. 
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\Aaj^{9i) . . . Aa„{9n)) are the multi-particle states relative to the excitations of mass m^, 
with relativistic dispersion relations given hj E = cosh 6, p = mk sinh 9, where 9 is 
the rapidity variable. The spectral representation (|4.36| ) is obviously an expansion in the 



parameter e « , where ^~ = mi is the lowest mass-gap. 

Basic quantities of the large distance approach are the Form Factors (FF), i.e. the 
matrix elements of the operators (pi on the physical asymptotic states 

FZ...aS0i,...,9^) = {Q\vmAaAOi),---,AaSOn)) • (4.37) 

It is worth emphasizing that the above quantities are unaffected by renormalization ef- 
fects since physical excitations are employed in their definitions. For scalar operators, 
relativistic invariance requires that the FF only depend on the rapidity differences 9i — 9j. 
Postponing a more detailed analysis of the analytic properties of the FF, let us first dis- 
cuss the behavior of the series (|4.36|) for the purpose of evaluating the integral hiR) in 
eq. 



First of all, it is convenient to order the multi-particle states entering the sum ( [4.36| ) 
according to the increasing values of the total sum of their masses E!^^''"'""" = Yl,k=i so 
that, for r ^ ^, the functions 5'„(r) behave as the decreasing sequence gnir) ~ e"**^"^ 
Apart from = (0|v9j(0)|0)(0|(^j(0)|0) relative to the VEV of the fields (which however 
does not enter the connected correlator), the first approximation of the correlator is given 
by the first term of the expansion ( |4.36| ) 

/ + 00 rlf) 
-oo Zn 



g -mircoshS 



= Ff ^ - Ko(mir) , (4.38) 

TT 

where -^'0(2;) is the modified Bessel function. Sometimes it may occur that the one- 
particle FF of the fields are zero for symmetrical reasons and, in this case, the leading 
approximation of the connected correlator is given by the function g2{f'), expressed in 
terms of the two-particle Form Factors Fl^^\^{9i — 6*2). Since (|4.36| ) is an exact expansion 

r 

in e ? , its truncated series to the lowest terms is expected to provide an accurate ap- 
proximation of the correlators in the interval r ^ C,- However, the convergence property 
of the truncated series is much better and as a matter of fact it neatly approximates 
the correlator up to the region r ~ ^, as has been checked in several examples (see, for 
instance |]10|, [Tl|, |I4|, ^). Therefore the truncated spectral series is assumed to 



estimate the integral l2{R) in eq. ( [4.24 ) with a reasonable confidence: obviously, the more 



terms included in the series ( [4.36| ) results in a better evaluation of the integral I2{R)- 
The problem is then to determine how efficiently we can assess the matrix elements of the 
order parameters on the asymptotic states. Let us discuss separately the cases when the 
off-critical theory corresponds to a non-integrable QFT or to an integrable one. 
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For a non-integrable QFT (as, for instance, the QFT resulting from the magnetic 
deformation of the TIM), unfortunately it is difficult to go beyond the one-particle Form 
Factors of the lowest particle states, {0\ipi{0)\Ak) . In fact, due to creation and annihilation 
events in the scattering processes of these theories, i.e. to the non-elastic nature of its 



S-matrix, the FF satisfy the infinite coupled set of Watson's equations 



F.^{9^), . . . , e„) = (0|yp(0)|A(^i) . . . A(0„))in = (4.39) 

CO „ jnl JQi 



m=0 



/■ fif)' fif)' 



m=0 



obtained by inserting the unitary sum on the out-states in the definition of the original 
Form Factor. Consequently, higher-particle FF have a non-trivial analytic behavior - 
with branch cuts at all production thresholds - which in practice precludes their exact 
determination. Hence, for non-integrable theories the best we can do is to estimate the 
large-distance expansion of the correlators only in terms of the lowest one-particle states. 
Moreover, the one-particle FF of these theories (0|(/3j(0)|Afc) cannot be determined by 
first principles and for their evaluation we have to rely on some numerical determinations, 
as discussed in next section. Although this situation may appear disappointing from a 
theoretical point of view, it is worth stressing that for all practical purposes, one can reach 
a reasonable estimate of the integral I2{R) also in the non-integrable case. This can be 
checked, for instance, by comparing the values of the integrals - obtained by approximat- 
ing the correlators by the lowest Form Factors - with their exact values obtained by the 
A-theorem, when the latter applies. 

For an integrable QFT, the situation is much better since, in principle, there is the 
possibility of determining exactly all Form Factors of the theory. For a detailed discussion 
of the calculation of the Form Factors in an integrable QFT, we refer the reader to the 



original literature [|T2, ^ ^ 0- Here we simply recall the basic equations of the two- 
particle FF |T(]| since, based on the fast convergence property expected for the spectral 
series [§, they will be the only terms employed in the following for approximating the 
correlators in their large-distance expansions (in addition, of course, to the one-particle 
ones). 

Assume, for simplicity, that the spectrum of the integrable QFT is made of the scalar 
particles Ai. Let Sab{0) be the elastic scattering matrix of the particles Aa and A},. In 
this case, the two-particle Form Factor F^f^{9) is a meromorphic function of the rapidity 
difference 6 satisfying the equations 

KM = Sa,{e) F^,{-e) , (4.4o) 
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F^,{zn + 6) = F^,izn - 9) . (4.41) 
Let F^^"'{9) be a solution of eqs. ( [4. 401) and ( |4.41| ), free of poles and zeros in the strip 



lm6 E (0,7r). By requiring asymptotic power limitation in momenta, F^fj{6) must be 
equal to F^™(^) times a rational function of cosh 6*, with the poles thereof fixed by the 
singularity structure of the scattering amplitude Sab{d). A simple pole in F^i^{9), like the 
one in Figure 1 induced by the simple pole of Sab{d) with a positive residue, gives rise to 
the equation 

F,^,(^^m:,)^-^%F,^ , (4.42) 

where is the on-shell three-particle coupling, also determined by the S-matrix. A 
more detailed analysis is, in general, required when the S'-matrix presents higher order 
poles (see the discussion in and Appendix D). The FF may also have kinematical 
poles which however do not appear at the two-particle level if the operator ip{x) is local 
with respect to the fields which create the particles. In conclusion, the two-particle FF 



can be expressed as [10 



Fm = #^i^r"(^) , (4.43) 

where Dab{d) and Qafe(^) polynomials in cosh^: the former is fixed by the singularity 
structure of Sab{0) while the latter depends on the operator V5(x). An upper bound on 
the order of the polynomial Q'^bi^) given in terms of the conformal dimension Ai^. In 
fact 

lim F'^(^) ~ e^^l^l , (4.44) 



with TO 



Z/^ < . (4.45) 

Further equations on the polynomial Q^bi.^) can be obtained when the field v^(x) is pro- 
portional to the trace of the stress-energy tensor T^u{x). In fact, as a consequence of 
the conservation law d^jT^^ = 0, the FF of Q{x) for two different particles Aa and Ab 
must contain a term proportional to the Mandelstam variable s = {pa + Pb)'^ of this state 
[ p!0| , ^ , so that it can be factorized as 

Moreover, in this case we have the normalization conditions for which reads 

FZ{^rc) = (AMlQmMOa)) = 2vrm^ . (4.47) 

The above discussion relative to the FF (and their generalization in the case of kink 
excitations) will be useful in Section 6 for the determination of the generalized suscepti- 
bilities for the integrable deformations of the TIM. 
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5 Numerical Methods: Truncated Conformal Space 
Approach 

The Truncated Conformal Space Approach (TCSA) has been introduced by Yurov and 



Zamolodchikov ||38| for a numerical evaluation of the non-perturbative effects relative to 
the off-critical models. It consists in studying the numerical spectrum of the off-critical 
Hamiltonian on a infinite cylinder of circumference i?, acting on the Hilbert space of the 
conformal states. Once a truncation at a suitable number of states is made, the problem 
reduces to perform a numerical diagonalization of a finite dimensional Hamiltonian. 
For the off-critical Hamiltonian we have 



H = Ho + V 



(5.1) 



where Hq is the Hamiltonian of the conformal fixed point on the cylinder and V = 

R 

Qij dvipi{w) where (fi{x) is one of the relevant perturbation (with w = u + iv is the 



coordinate along the cylinder and the tilde indicates quantities defined on the cylinder). 
By using the conformal transformation z = e"R the conformal theory on the cylinder is 
mapped onto a plane and therefore Hq can be expressed in terms of the usual conformal 
generators Lq, Lq and the central charge c [pT| : 

2n 



-(Lo + Lo-- 



27r 



R 



2A,: 



(5.2) 
(5.3) 



The spectrum of H depends on the dimensionless parameter QiR?'"^^^ and the value 
of Qi can be fixed^ such that the mass gap is equal to 1. Let be the inverse of the 
matrix hmi = {m\T) introduced to account for the nonorthogonality of the conformal basis. 
Denoting by 
have 



h^^li (J.\H\n) the matrix elements of the perturbed Hamiltonian, we 



Hr, 



27r 

'r 



R \ 

(2A„ - C/I2)5mn + 27r(7, ( — j 



2{1-Ai 



(5.4) 



The matrix elements {l\ipi\n) can be computed in terms of the structure constants of the 
OPE and the action of the conformal generators L„ on the states. For the perturbation 
of the minimal models of CFT a numerical algorithm has been designed to compute the 
above matrix elements and to perform the diagonalization of the off-critical Hamiltonian 
by including the conformal states and their descendants up to the fifth level of the Verma 



module |52|. In the case of the TIM this is equivalent to truncate the number of states 



By using this normalization, we were able to determine, in particular, the constant Ci entering the 



relation between the coupling and the mass in the magnetic tpi deformation (see eq. (4.4)) 
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to 228. Once the Hamiltonian H has been diagonahzed for different values of R, one 
can extract the spectrum as a function of R, in particular the low energy eigenvalues (see, 
for instance Figure 2) and also their associated eigenvectors. In this way, it is possible to 
determine the masses of the lowest particles, several vacuum expectation values, some of 
the one-particle Form Factors and also some of the generalized susceptibilities. There are 
however certain limitations of the method, one of them due to the truncation performed 
in the number of the conformal states employed in the algorithm (see the discussion 
p7|). In fact, all the quantities of interest are infrared data, i.e. relative to the 



m 



dynamics of the system on the cylinder in the limit i? — > oo, which is however dominated 
by truncation effects. This means that in order to extract reliable infrared data one has 
to look at the spectrum within an interval of R sufficiently large but still unaffected by 
truncation errors. This interval will be called the "physical window" . Another limitation 
of the TCSA occurs when the conformal dimension A of the perturbing operator is such 
that A > i. In this case, in fact, the renormalization of the operator prevents to reach a 
proper scaling behavior of the energy levels and the only quantities which can be extracted 
with reasonable confidence are the energy differences AEn{R) = En{R) — Eo{R). 

5.1 Vacuum Expectation Values by TCSA 



As shown in ref . , the knowledge of the eigenvectors in the Truncated Conformal Space 
Approach allows a numerical estimation of the Vacuum Expectation Values of several 
order parameters. These quantities are defined by the limit 

Bji= lim {6\0j\6) i\gi\-^^ , (5.5) 

it— »oo 

where |0)j is the vacuum (on the cylinder) relative to the off-critical theory along the i-th 
deformation. On the other hand 

(01^.10). = [-^) (0|c,,|0). = (^-j ^c,^_^o ' (5-6) 

where ip^ designs the m^^ component of the ground state vector expressed in terms of 
the conformal basis. Due to the truncation effects discussed above, the limit i? — > cxd in 
eq. (|5.5| ) in practice means that one has to consider R large enough such that the VEV 
reaches a saturation plateau. This saturation can be numerically controlled by requiring 
that {0\^j\0)i ~ R'^^', I.e. 

1 rfln(0|(^,|0), 



2Aj d\n R 



1. (5.7) 



By using this procedure, several VEV for different deformations of the TIM were deter- 
mined in We have reproduced and confirmed the results of |^ for the two most 
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relevant perturbations a and e of the model (see Table 0). In Table [1^, we have also 
included the exact VEV extracted from ref. of the TIM perturbed by the thermal 



operator (for both high and low temperature phases) in order to test the feasibility of 
the numerical approach. As evident from this table, lower the dimension of the operator, 
better the accuracy of the method. This is easy to understand since the computation 
of the VEV on the cylinder is equivalent to compute the VEV at a finite-temperature, 
a situation analyzed in [^. The numerical determinations of the VEV's relative to the 
</93 and (/94 deformations turn out to be quite inaccurate for the renormalization reasons 
discussed above (A3 ~ \ and A4 > but their exact values can be nevertheless extracted 
from the results of ref. 



5.2 Numerical Determination of the One— Particle Form Factors 

As discussed in Section 4, once the Form Factors of an operator are known, its correlation 
functions can be written as an infinite series over multi-particle states, eqs. ([4.35|) , (|4.36| ), 
and the restriction of these series to the first one-particle states already provides a reason- 
able approximation of their long distance behavior. These one-particle matrix elements 
(0 I V5i(0) I y4fc) can be numerically determined along the lines followed for estimating the 
VEV of the various operators. Namely, one has to replace the infinite volume vacuum 
state (0 I with its components relative to the eigenvector of the truncated Hamiltonian, 
analogously for the vector relative to the one-particle state, and then use the matrix 
elements of the field yjj in the truncated basis. The resulting quantity finally needs to 
be multiplied by y/mkR because the normalization of the one-particle states on a finite 
volume differs precisely for this factor from the one in the infinite volume. The quantity so 
determined, plotted versus the radius R of the cylinder, presents in many cases a plateau 
in the region of the physical window which therefore provides its numerical estimation. As 
an example of such determination see Figure 3 where the matrix element (0 | V9i(0) | Ai)i 
is plotted as a function of R in the case of the leading magnetic deformation of the TIM. 
A plateau is clearly reached for i? ^ 12. We have performed this numerical calculation 
for all excited states under the two-particle threshold relative to the lowest mass gap 
for the first three deformations. The results are in Tables 0, ITTl O and ITTI, in units of 



the opportune power of the coupling constants. As for the numerical determination of 
the VEV's, also in this case the lower the dimension of the operator, the accuracy of 
the method improves. Moreover, the numerical errors are usually larger for those matrix 
elements involving states which are closest to the threshold. 

Let us conclude this section with a general remark on the one-particle Form Factors. 
We have already discussed that for the non-integrable deformations the knowledge of these 
matrix elements is crucial in obtaining at least a non-trivial estimate of the correlators 
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and their determination necessarily passes through a numerical approach. However, even 
in the more favorable case of integrable theories, it may occur that the determination of 
the one-particle FF can be only obtained by a numerical approach. Despite the existence 
of a manageable set of recursive equations which link the various n-particle Form Factors 
in the integrable models, the solutions of these recursive equations need an initial input 
which cannot be often obtained even by employing the cluster property of the Form Factors 



Q Under this circumstance one has to necessarily resort to other methods for 
obtaining the one-particle FF's and the TCSA may help in this respect. An explicit 
example of the situation discussed above is provided for instance by the Form Factors of 
the operators ipi in the thermal deformation, as discussed in the Appendix |^. 

5.3 Numerical Computation of the Susceptibilities 

The TCSA also allows a direct numerical estimation of the susceptibilities. Since they 
are defined as 



2 £ '^1+'^ 



= r}fc^7, , (5.8) 

one needs to numerically evaluate the derivatives of the VEV's with respect to the different 

R 

couplings. Hence, a small perturbation g^^ J dv(pk{w) is initially added to the Hamiltonian 



( |5.1| ), with the values of the coupling constant gu chosen in such a way to alter the 
spectrum of the unperturbed theory only of small percent. To express the final result in 
unit of \gi\ ^-^^ , it is convenient to write the coupling constant gt as 

gk = akigl'^' . (5.9) 

The next step consists in computing the expectation value of {0j)i^k as in eq. (|5.6|) by 
varying aki- As our typical sampling, we have considered 5 — 10 different values of aki and 
then we have extracted the numerical estimates of various susceptibilities by a linear fit 
of the data in the physical window of the R axis where the VEV presents a plateau. The 
data relative to the two strongest relevant deformations are in Tables |T^ and |15[ Their 
values are reasonably close to the ones obtained by the fluctuation-dissipation theorem or 
to their exact values, when available from the A-theorem sum rule. The only exceptions 
are those relative to the susceptibilities relative to the ip^ operator, where there is a 10% 
mismatch. A non trivial and internal check of our estimates is given by the symmetrical 
relation F*^ = F^^- shown by the data. As for similar calculations discussed above, this 
method seems however to be inappropriate for the higher dimension perturbations relative 
to the fields ips and ip^. 



24 



6 The Four Relevant Perturbations of the TIM 



In this section, we will discuss in some details the field theories relative to each individual 
relevant deformation of the TIM, i.e. those associated to the actions 

Af = AcFT ±gijd?x^i{x) , z = 1, . . . , 4 (6.1) 

The qualitative form of the effective potential relative to the different off-critical defor- 
mations of the TIM is shown in figure 4. In addition to the spectrum of the off-critical 
excitations, for each field theory we will present the main formulae involved in our esti- 
mation of the universal amplitude ratios of this model. Some aspects of the calculations 
of the correlation functions and the VEV are also discussed, refering to the appendices 
for all technical details. Here it is worth to comment on an interesting feature of the 
Conformal Perturbation Theory common to all deformations, i.e. the appearance of log- 
arithmic terms both in some VEV and in the calculation of some susceptibilities. The 
origin of some of these term can be traced back to the existence of some peculiar resonance 
conditions involving the anomalous dimensions of this model. The first of them is given 
by 

Ai = A3 + A4 - 1 , (6.2) 

which is equivalent to say that the scale dimensionality of the coupling gi equals the 
product of (73 and 5(4, namely gi ~ g^g^- It is easy to see that this resonance condition may 
infiuence the calculation of some susceptibilities. In fact, even though each individual off- 
critical field theory is defined by the one-coupling action ( |6.1| ), nevertheless the calculation 
of the susceptibilities r*^ requires to consider initially the multi-coupling action 



A = Af + gj J (fx (pj{x) + gk J d'^x ipk{x) , 



(6.3) 



and to take the limit gj = gk = only at the end of the calculation. Hence, the resonance 
condition gi ~ (73(74 may spoil the naive form of the above action ( |6.3|) with the presence 
of additional terms. Explicit examples of this phenomenon are commented in the next 
sections. Another resonance condition which also infiuence some of the calculations is 
given by 

Ae = 1 + A4 - A2 , (6.4) 
where Ag = | is the anomalous dimension of the irrelevant field e" . 

6.1 The Magnetic cpi Deformation 

This is the most relevant deformation of the TIM and the only non-integrable one. Hence, 
most of the results relative to this deformation are obtained by the help of the numerical 
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approach. The numerical analysis of the spectrum, first performed in shows that 
there are two different one-particle states with mass ratio 1712/ rrii ~ 1.61 (see Figure 2). 
By setting the value of the first mass to be 1, its relationship with the coupling constant 
gi is numerically determined to be 



40 40 



mi = Cigr ^ 3.242... . (6.5) 

The VEV of the different fields have been numerically computed and their values are in 
Table |^. By applying eq. ( 4.28 ), it is easy to check that there is no UV mixing of the 



operators for this deformation and therefore no need for their UV renormalization. In 
order to compute the various susceptibilities, we have decomposed / dPx {ipi{x)(pj{0))i 
into the two integrals /i(-R) and h^R), as discussed in Section 4. The UV part of the 
correlator has been approximated by the short range expansion ( [4.34|) with the employed 
values of {Cfj)i^^ reported in Table |^. They were computed as explained in appendix B. 
For the IR part of the correlator, the non-integrable nature of this deformation forces us 
to truncate the spectral expansion to the one-particle contributions only 

{^.{x)vM) ^ j:FrFpK,{m,\x\) . (6.6) 
1=1 

The numerical estimation of the one-particle Form Factors Fi, expressed in opportune 
units of gi, can be found in Table |. 

As a concrete example of the above procedure, let us consider the correlator {ip2{x)ip2{0))i- 
Its UV expansion reads: 



-4A 



+ (CJ.)S"B,. (^|^)% 0((™,r)V (B.)^ (^) 



4A2 



(6.7) 



The above expression is expected to provide an accurate approximation of the correlator 
up to rriir ~ 1. Indeed, in a plot of the UV and the IR approximations of this correlator 
(Figure 5), a satisfactory overlap between the two curves is observed around mir ~ 1: 
this makes us confident on the estimation of the susceptibility extracted by integrating 
the above correlator. The same situation occurs for the other susceptibilities and the final 
results are in good agreement with those extracted by the A-theorem sum rule or their 
direct numerical estimation by TCSA. All these data are reported in Table ||. The only 
exceptions consist in the calculation of the susceptibilities and for which some care 
is required due to the anomalous dimensions of the fields involved and to some subtleties 
in the conformal perturbation expansion of Cfj{g;x). In these case, for instance, the 
naive numerical integration of the corresponding correlators cannot be performed because 
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of their divergences at r ^ 0. Let us first consider r44. A natural way to regularise the 
integral / (Px {Lp4(x)ip4{0))i which near r — goes as 27r / drr'"^^^, consists in introducing 
a UV cut-off a so that its UV divergence is easily extracted 

= — + finite part . (6.8) 

as 

Once this divergence is subtracted, the finite part of the integral, being cut-off indepen- 
dent, may be regarded as the actual regularized susceptibility. For the magnetic defor- 
mation its value is in Table ^ with 5 — 10% of approximation and this quantity can be 
used later on in the evaluation of the universal ratios. The same strategy has been also 
adopted for the evaluation of the susceptibilities r44 relative to the other deformations of 
the TIM since the divergence of the integral is simply due to the conformal properties of 
the field fii^x) and does not depend on the particular deformation considered. 

The situation is different for r34. In this case, by using the operator product expansion, 
it is easy to see that the integral / d'^x {(p3{x)ip4^{0))i behaves in the UV region as 

. (6.9) 

Ja r 

and therefore it presents a logarithmic divergence 

where G is the finite part. When Bn is different from 0, G is not uniquely defined since 
it varies by changing the value of g^. Hence, contrary to the previous case, the finite part 
of the integral cannot be used to define universal ratios although the amplitude in front 
of the logarithmic term is an unambiguous quantity in QFT which may enter universal 
combinations. It is worth pointing out that this situation is not peculiar of the TIM but it 
is already familiar in the context of studying the specific heat dependence in the standard 
two-dimensional Ising model (see, for instance |12|). Obviously the above considerations 
apply for all the susceptibilities r34 relative to the other deformations of the TIM. 



6.2 The Thermal (p2 Deformation 



This is an integrable deformation of the TIM |21|, |22[. When g2 > 0, the coupling 
to the thermal field moves the TIM in its high-temperature phase, where a unique Z2 
symmetric vacuum state is present. When g2 < 0, we reach the low-temperature phase 
of the model, where there is a spontaneously breaking of the Z2 spin symmetry and 
therefore two degenerate symmetric vacua. In the high-temperature phase there are 
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ordinary massive particle excitations whereas in the low-temperature phase there are kink 
excitations and bound states thereof. The two phases are related each to the other by a 
duality transformation. The off-critical model possesses higher conserved charges whose 
spins are s = 1, 5, 7, 9, 11, 13, 17 (module 18), i.e. the Coxeter exponents of the exceptional 
algebra Ej. The integrable structure of this deformation originates from the conformal 
decomposition (p.l3|) together with the pairing of the energy operator e(x) with the adjoint 
representation of the WZW model on (£'7)2 [^. The existence of an infinite number of 
conservation laws implies the elasticity and factorization of the scattering amplitudes. 



These amplitudes have been computed in ^T], and their concise expressions can be 



found in Table 2 of ref. The exact mass spectrum of the excitations can be extracted 
from the pole structure of the S matrices (see Table |^). With respect to the Z2 spin 
symmetry of the model, in the high-temperature phase there are three Z2 odd particle 
states (the ones relative to the masses mi, and mg) and four Z2 even (those relative to 
the masses m2, m.4, and my). In the low-temperature phase, the three Z2 odd particle 
states become kink excitations interpolating between the two degenerate ground states 
whereas the four Z2 even ones play the role of breather states. This leads, in particular, 
to an interesting prediction on the universal ratio of the correlation lengths above and 
below the critical temperature. In fact, if the correlation length is defined according to 
the leading exponential falling off of the spin-spin connected correlation function in the 
limit I X I 



X 



(0|¥.i(x)y.i(0)|0),±~exp(^-^j , (6.11) 

(where the indices ± refer to the high and low temperature phases respectively), from 
the Z2 symmetry property of the a field, the self-duality of the model and the spectral 
representation of the above correlator we have 

= = 2 cos — = 1.28557... (6.12) 

i- mi 18 ^ ' 

For this deformation the relationship between the mass-gap mi and the coupling 

5 

constant §2 is given by mi = €292 where 




47r2r(^)r(^)n'^'' 

C2 = ( 1 I % ^ = 3.7453728362 . . . (6.13) 

The VEV's of all relevant operators, both in the high or in the low temperature phases, 
have been computed in [5D| and their values successfully compared with their numerical 
determination the expressions of the expectation values of the operators a and e 
are finite, whereas those of the operators a' and e' are naively divergent (see eq. (^3)) 
and need therefore a regularization (see Appendix C). They are presented in Table |10|, 
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in opportune units of the associate power of the couphng constant. Concerning the UV 
properties of the fields, eq. ( [4.28| ) predicts that there is no UV mixing of the operators for 
this deformation and therefore no need to implement their UV renormalization. 

Let us now turn the attention to the computation of the various susceptibilities by 
following the strategy explained in Section 4. The values of the coefficients {Cfj)2^^ en- 
tering the UV expansion of the correlators are in Table 0. As a concrete example of 
these calculations, let us discuss the susceptibility amplitude F^^ relative to the correlator 
{ipi{x)ipi{0))2- Its UV expansion reads: 



(V5iW^i(0))2 = r 



-4Ai 



6 , , 9 



V L2 / V ^2 / V '-'2 

\ 4Ai 



3\ 



(6.14) 



The IR part of the correlator is approximated by taking initially into account the first 
four one-particle states of the spectral expansion of the correlator (see eq. ( |4.36| )) 

{Vi{x)^m ^j2FrFpKo{mi\x\) . (6.15) 
1=1 



The numerical values of the needed one-particle Form Factors are in Tables |Tl| and 112 



relative to the high and low temperature phases respectively. Although these matrix 
elements are already able to reproduce with a reasonable accuracy the infrared part of 
the correlator, due to integrability of this deformation some two-particle Form Factors are 
also available and therefore they can be used to improve the estimation of the correlators. 
Their calculations, together with some subtleties which occur in this case, are discussed 
in Appendix D. The Form Factors of the operator e(x), which plays the role of the trace 



of the stress-energy tensor for this deformation, were computed in ref. [47 



The above strategy has been applied for the estimation of all the correlators. An over- 
lap between the UV and IR approximations of the correlators has been usually observed 
in the region rriir ~ 1, which may be regarded as a consistent check of our approach. 
Such overlap is shown in Figure 6 for the correlator {ipi{x)ipi{0))2- 



In closing this subsection, some comments are in order for Table |T3| which collects 
the values of the first corrections to the structure constants in the thermal deformation. 
Three of these coefficients contain some logarithmic dependence which however do not 
particularly infiuence the numerical approximation of the short-distance of the correlators 
since these are higher-order corrections. The presence of these logarithms are due to a 
resonance phenomenon between the conformal families, i.e. it occurs when the conformal 
dimensions of two operators differ for an integer number of times 1 — A, where A is the 
dimension of the perturbing operator. This situation is encountered here because 

A6-A4 = (I-A2). 
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Hence, we have for instance {Cf-^)2 {mr) = ai(ln(mr) + bi{m)) with the coefficient ai 
that can be computed by using the methods explained in Appendix 2. Since its algebraic 
expression is rather cumbersome, we prefer to report here just its numerical value ai = 
—0.1510653.... While ai is an unambigous quantity in QFT, the other term bi{m), 
on the contrary, is a scale dependent quantity. Similar resonance phenomena are also 
encountered in the sub-leading magnetic deformation of the TIM, which is the subject of 
the next subsection, where they play a more important role since they contribute to the 
lower orders in the UV expansion. 



6.3 The Sub— leading Magnetic (f-^ Deformation 

This is an integrable perturbation [^, It is generated by the less relevant magnetic 
field and obviously breaks the Z2 spin symmetry of the critical point since the field <^3 is 
a Z2 odd operator. The resulting massive theory presents some interesting features, as 
ffist outlined in ref. |2^. The theory presents two degenerate albeit asymmetrical vacua 



(denoted by IO2) and IO4)). There are two massive kink excitations of mass m and one 
breather bound state thereof with the same mass. The exact S'-matrix of this theory has 
been computed and analyzed in The relationship between the coupling constant 

8 

and the mass gap of the theory is provided by m = Csg^ with the constant C3 given by 



37] 



V3r(i)r(|; 
-r(|) 



4.92779064... (6.16) 



The VEV's of relevant operators (pi, ip2, V^s have been exactly computed in The 
expression of the VEV's (0s|v34|0s), s = 2,4 given by eq. ( [4.9|) are a priori divergent 
and need to be regularized. As shown in Appendix C and also discussed below, they 
present a logarithmic dependence on the coupling constant g^. The complete set of VEV, 



expressed in unit of the appropriate power of the coupling constant, are in Table |T6 
Notice that, according to the formula ( |4.28|) , the only field which requires an ultraviolet 
renormalization is precisely the field (p4,{x) which mixes (logarithmically) with the field 
ipi{x) under perturbation theory. 

The application of the A-theorem sum rule allows the exact determination of the four 
susceptibilities T^j (as already discussed in Section 6.1, the susceptibility T^^^ contains 
a logarithmic dependence on the coupling constant g^). Also for this deformation we 
have followed the general strategy explained in Section 4, i.e. we have ffist pursued 
the matching between the UV and the IR expansions of the correlators and then we have 
performed their integration for extracting the susceptibilities. Let us discuss the correlator 
{ipi{r)ipi{0))3 to exemplify some new difficulties arising in the perturbative evaluation. Its 
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UV expansion is given by 

/ 1 6 

/ N / N> A A fmr\5 fmr\5 / mr 



+ (Cmr)B,,{^)\...-iB,,r{^)"') . (6.17) 

There are two different sources of problems in this expansion. First of all, the VEV {ip^)^ 
contains a logarithm dependence on the coupling and a UV regularization is needed. By 
using eq. ( |6.10| ) and the A-theorem sum rule 

B,s = ~^-^^Tl, (6.18) 



it can be easily shown that 

where e is some arbitrary UV cut-off and A ~ ,^3. Secondly, the coefficient (Cj^i)3^^ 
contains some logarithmic divergences due to resonance problems already encountered in 
the thermal perturbation and cannot be calculated by using the Mellin transformation 
method in the usual manner. However, the Mellin transformation method can be properly 
generalized to obtain the correct logarithmic term, as we show in the following. {Cl^)'^\r) 
is defined as usual by the regularized integral (notice that we have restored the dependence 
on r) 

{C'uYi\r) = - f d'z{v,{oo)v,{r)^,{z)^,{0)) ■ (6-20) 



Using notation of Appendix A, {Cli)'^\r) may be regarded as the regularized limit s ^ 
of 



G{2 -s;r) = / d'z \z\'' {Moo)MrW3iz)M0)) 

d'zd\\u\^''\u - l\''\u - z\^''\z\^''-'\z - (6.21) 



where the coefficients a, 6, c, d, e are expressed in terms of the weights by the Coulomb 
gas formalism and s is the Mellin transformation parameter. After putting r = 1 in the 
previous equation, we have 

— Ci 

G(2-s;l) = -^ + 0(l) , (6.22) 

s 

i.e. there is a first order pole which does not allow us to calculate (Cj^i)^^'* in the usual 
manner. However, by multiplying the previous expression with the expansion (see (|C.48|) ) 

r-' = m\mr)-' = [I - slnmr + 0{s^)\ , (6.23) 
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where m is an arbitrary mass scale, we obtain the expansion in s 

G(2-s;r) = -^ + — cilnmr + 0(l) . (6.24) 

s 8 

Hence, (C*ii)3^'*('") may be taken as the regularized version of the previous expression, i.e. 
the one obtained by discarding the simple pole divergence 

(Cn)^^^(r) = ^ci(ln(mr) + const) , (6.25) 
8 

where the constant term clearly depends on the scale m. By looking at eq. ( |6.21| ) and by 
taking into account the dependence of a, b, c, d, e on the weights of the primary fields, it is 
easy to see that a shift of the parameter ^, defined in eq. ( |C.31|) , ^ = -rz ^ + ^ would 

6 

have led to the same result, with e playing the role of s. The two terms in ^ lead in 
the correlator to ^ 

yci^ia (^) ' (In(mr) + A{m)) , (6.26) 

where A is a cut-off dependent quantity which can be determined by requiring the match- 
ing between the infrared expansion of the correlator with its UV part. The infrared part 
is approximated by the one-particle Form Factors of the magnetic and thermal opera- 
tors (which are reported in Table 0). There is a non-trivial consistency check of the 
above procedure. Indeed, the correlation function can be computed in the two differ- 
ent vacua IO2) and IO4). Once the quantity A has been fixed by computing the correlator 
(02|v^i(^)v^i(0)|02), the same A should also work for the other correlator (04|v9i(r)(^i(0)|04), 
as it is indeed the case. 

The above procedure has been also employed for the correlators {^P2{r)^P2{0))3 and for 
{ipi{r)ip2{0))3. For this deformation, however, we were unable to reach any definite result 
on those correlators involving the ip^ operator because its one-particle FF were completely 
inaccurate, thus preventing a reliable evaluation of the infrared part of the correlators. 
The results for all the susceptibilities we were able to compute for this deformation are 



in Table M 



6.4 The Vacancy Density (^4 Deformation 



The field theories originated by the perturbation t{x) = ip4{x) with (74 > and (74 < 
are both integrable [^. The most elegant way to get an insight on these quantum field 

i|, E6l. Since the 



theories is to use the supersymmetric formulation of the TIM [|T9 



field t{x) is the top component of the superfield ( |2.10| ), the off-critical dynamics may be 
described by the action 

■1 



A 



-VAfVAf + Af^ + g^Af 



(6.27) 
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After eliminating t{x) by its algebraic equation of motion, the interaction terms of the 
above lagrangian are given by 

i^^e+l(le^ + gX • (6.28) 



2 V2 

Hence, for > the ground state energy is nonzero and supersymmetry is spontaneously 
broken: the scalar field acquires a mass whereas the fermionic field remains massless and 
plays the role of goldstino. This is nothing but the Majorana fermionic field of the fa- 
miliar two-dimensional Ising model which is in fact the ending point of the massless 
Renormalization Group flow originating from the g4 > deformation [^. On the mass- 
less Majorana fermion of the Ising model supersymmetry is implemented non-linearly. 



The exact massless ^-matrix has been computed in |]2^ and the related massless Form 
Factors in p5[. The crossover phenomena, where the exponents characterizing the leading 
singularity of an observable change from its tricritical to its critical value, have been stud- 
ied in Along this flow, the conformal dimension of the magnetization field changes 



from its tricritical value to ^ of the Ising fixed point, the conformal dimension of the 
energy operator varies from jq to j and finally the conformal dimension of the vacancy 
density operator changes from ^ to 2, since this operator becomes - at the end of the 
Renormalization Group flow - a descendent of the Identity operator. The sub-leading 
magnetization operator also becomes at the end of the RG flow a descendent of the mag- 
netization field of the Ising model. For the massless nature of this theory we will not 
compute the relative universal amplitude ratios. 

For g4 < the ground state energy vanishes therefore supersymmetry is unbroken and 
both the scalar and the fermion fields become massive. As argued in [^, to describe the 
dynamics of the system it is more convenient to adopt the usual Landau-Ginzburg poten- 
tial in terms of the order parameter a{x). This potential presents three-fold degenerate 
vacua, labelled by |0_i), |Oo), |0+i), where |Oo) corresponds to the disordered vacuum and 
|0±i) to the two disordered vacua, symmetrically placed with respect to the origin (see 
Figure 4). The elementary excitations in this phase are massive kinks, which interpolate 
between the ordered vacuum and the disordered vacua and are denoted by Ko±i or K±io. 



The associated scattering theory has been discussed in |26, 56 



Let us consider the amplitude ratios for this deformation. The VEV of {ipi)4 and {1^2)4 
can be directly computed by using eq. ( [4.11| ), since the integral converges. The VEV's of 
!f3 contains a logarithm divergence and is calculated in Appendix C. This is the mirror 
situation encountered in Section |6.3| for the VEV's (^54)3 and therefore it can be fixed as 
in the previous subsection: 

(</^3)4 = ^5i4 (^In ^ + const^ • ^^-^^^ 
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Finally, by supersymmetry we expect {^i)^ = 0, because the field fi^x) plays the role of 
the trace of the stress-energy tensor for this deformation and indeed this is in agreement 
with formula ( [4.11| ). The values of the VEV may be found in Table 



Concerning the UV properties of the theory, there are two fields which need renor- 
malization. The first is the sub-leading magnetization (t'{x), which mixes logarithmically 
with the magnetization field a{x). The second is the vacancy density field t{x), which 
mixes with the energy density e{x). To avoid the difficulties related to the mixing in- 
duced by renormalization, for this deformation it is convenient to rely only on the Form 
Factor expansion for estimating the correlation functions. In fact, all fields, except the 
sub-leading magnetization turn out to be uniquely identified by their symmetry 

properties. Moreover, the spectral series based on the Form Factor are also able in this 
case to capture successfully the ultraviolet behavior of the correlators Hence, for this 
deformation the integrals ( |4.24| ) will be estimated only in terms of /2(0), i.e. / -^2(0). 
The two-particle matrix elements on the kink states for the operators ipi, and have 
been calculated in ||5^ by using a mapping of the TIM onto the dilute q-state Potts model, 
with q = 2. All these operators are coupled to states with zero topological charge, i.e. 
to those multi-kink states which begin and end to the the same vacuum. There are two 
types of such two-kink Form Factors, depending whether the vacuum^ is |Oo) or |0-|-i) 

F^U0i-02) = {0o\v^\Ko±l{el)K±^o{02)\0o) ; (6.30) 
^^10(^1-^2) = (0±i|<^,|i^±io(^i)i^o±i(^2)|0±i) , (6.31) 

where (pi designs the operator under consideration. The Form Factors can be conveniently 
parameterized as 

(6.32) 
(6.33) 
(6.34) 
(6.35) 

and their expressions are given by 

Fr (^) = -mm^e±i(-+^^) . '"t ^ • ^ ^o(^) , (6.36) 

smhj[U — ITT) 

4 smh ^{o — m) 

^..(,,.^^i^__^T,Wf„(.), (6.38) 





= F^^{9) , 


z = 2,4 




= Ff (e) , 


z = 2,4 




= ±F^'{e) , 




FlUe) 


= ±Fr(0) , 





^"Note that we use a different notation for the vacua compared to the one used in |4^. These notations 
enable us to specify more clearly the disorder vacuum |0o) (with zero magnetization) and the two ordered 
vacua |0±i) (with ±1 magnetization). 
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with 



2 JO 



a; sinh 2x cosh | sinh x 

2/ 



T (B) - P.n /9 r ^^^^(i - sin^(2m-g)^ ] 
T4^)-exp|2y^ -^hlhf sinh 2a; / ' 

T_(^^) = T+(^^ + 227r) , 
7 = — In 2 . 

With the knowledge of the first Form Factors, the spectral representations of the correla- 
tion functions involving the fields (pi{x), (p2{x) and ip4{x) are approximated by 

r lif) rif) 

(Oo|(/p.(x)y.,(0)|Oo)^ E / -7^^F^iA0i~d2)F^U6,-6{)e-\^\''^ ; (6.39) 

, J8i>92 2n Zn 



k=±- 



(o±i|vp.(x)^,(o)|o±i) ^ / -^^F^{o{ei-e,)Fli,{e,-e,)e-\^\^^ , (6.4o) 

where i,j = 1,2,4 and E2 = m(cosh 61 + cosh 62) is the energy of the two-kink asymptotic 
state. By integrating the above expressions, one can obtain the associated amplitudes. 
The convergence of the above spectral series has been successfully checked against the 
A-theorem sum rule (when this applies) . Our results for the amplitudes are in Table 
I2TI. 



7 Conclusions 

In this paper we have exploited some powerful techniques of Quantum Field Theory 
in order to compute an ample set of universal amplitude ratios for the scaling region 
of the two-dimensional Tricritical Ising Model. The determination of the thermody- 
namical amplitudes entering the universal ratios has been obtained by combining exact 
non-perturbative results coming from CFT (UV theory) and from Scattering Theory (IR 
theory). More specifically, we have used the ultraviolet data provided by CFT for setting 
up an Operator Product Expansion (and computing the first order approximation to the 
structure constants), and for finding eigenvalues and eigenvectors of the off-critical hamil- 
tonian by a numerical approach. In addition, we have used non-perturbative approaches 
related to the integrability of several deformations of the model to obtain important in- 
frared data, as for instance the exact values of the Vacuum Expectation Values of the 
order parameters, the exact spectra of the massive excitations and the first Form Fac- 
tors. A judicious use of the ultraviolet and infrared properties of the various fields which 
span the scaling region of the model have allowed us to reduce significantly the analytic 
efforts necessary to compute the universal ratios. Some of these quantities can be found 
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in Tables 22-26|^ As already pointed out in the text, this large set of quantities may 
be quite useful for an experimental investigation of the critical properties of this class 
of universality and we hope that the results presented in this paper may stimulate such 
experimental activity. It would be equally interesting to extend the theoretical approach 
discussed here to other two-dimensional models in order to reach a full control of the 
classes of universalities of low-dimensional systems. 
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In the calculation of the universal ratios, we have used the exact values of the susceptibilities 
when available from the A-theorem sum rule, otherwise we have used the arithmetic mean of their 
determinations obtained by the fluctuation-dissipation theorem and the numerical TCSA. 
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A Mellin regularization scheme 



The aim of this appendix is to discuss a regularization of the integrals ( |4.33| ) . They can 
be written as 

Cfj'^ = - f'd'wg{w,w) , (A.l) 



with 

g{w, w) = {AP{ooMw)^,{l)^,{0))cFT. (A.2) 

A regularized version of the divergent integrals (|4.33| ) is provided by the following function 
of the large distance cut-off R=l/a 

I{a) = — I d^w Q{a\w\) g{w,w) , (A. 3) 



where the cut-off function B(t) has a fast decreasing behavior at t — +00 to make the 
integral convergent and is equal to 1 for t ^ 0^. Of course /(a) diverges for a ^ 0"*" in 
a way that depends on the particular choice of but it converges to a finite value 

for a +00 thanks to the fast decreasing behavior of the cut-off function. The finite 
part of the integral ( |A.3| ) is independent of the parameter a and furnishes its natural 
regularization. It coincides with the analytic continuation of the integral ( |A.1| ) from those 
regions of the conformal weights for which it converges. The finite part of the integral can 
be explicitly calculated by means of the Mellin transform of the complex function J(a) 
defined by 

I{s) = / — a^J(a) . (A.4) 
JO a 

In fact, when the Mellin transform have simple poles, it provides the asymptotic expansion 
in powers of a of the function I{a) according to the formula 

I{a)=J2Res[a-'i{s)l=s, , (A.5) 

i 

where the sum runs over the poles. In our cases the above sum is finite since the theory 
presents only a finite number of IR divergent terms. The finite part of I{a) is therefore 
given by 

/o = Res[J(s)],=o = lims/(s) . (A.6) 

Finally, we also need the following theorem on convolution: if the function I{a) has the 
form of a convolution (as in equation ( |A.3|) ), then its Mellin transform is given by 

i{s) = -Q{s)G{2-s) , (A.7) 

where 6(s) is the Mellin transform of 0(t) and G{s) may be considered as the Mellin 
transform of the angular integral of g{\w\,arg{w)) with respect to the radial coordinate 



\w\ 



G{2-s) = d^w\w\-'g{w,w) . (A.8) 
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In our calculation we have used for 0(t) the function 

Q{t) = e"* , (A.9) 
whose Mellin transformation is the Gamma function 

Q{s) =r(s) . (A.IO) 
For the calculation of G{2 — s) we refer the reader to the next appendix. 

B Calculation of the {Cfj^^)k 

In this appendix, we show how to compute the first correction to the structure constants, 
i.e. the finite part of 

(CW = I (¥^p(oo)$,(^)v.,(l)y.,(0)) £w . (B.ll) 



The conformal four-point correlation functions entering the integral may be computed by 
means of the modified Coulomb Gas methods P]. In this approach the central charge c 
is parameterized by 



c=l-2Aal ; a± = ao ± y al + 1; (B.12) 
a+a- = —1 . 

The vertex operators are defined by Vnmi^) ='■ e*""™**^^-* : where $(a;) is a free scalar field 
and the charges anm defined by 

"nm = ^(1 -^)«- + ^(1 -"^)a+ • (B.13) 
The conformal dimension of the operator Vnm{x) is given by Anm = —oi^oinm with 

a;;;;;^ = 2ao - anm = ^(^ + ^)^- + ^(^ + • (B-14) 

The integrals encountered in the our computation were of two types. 

In the first case, no screening charge is needed, therefore they can be computed in a 
straightforward way. For instance, this is the case of the integral of the 4-point correla- 
tion function {ipi{xi){p2{x2)ip2{x3)ip3{x4)) . As an example, we provide the calculation of 
(^22)1^'' magnetic deformation of the TIM: 

(CsT^)! = - / (V^3(oo)y.2(l)¥^l(^)^2(0)) . 
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Since 2ai2 + 0:22 + «2i = 2ao, no screening charge is needed. Therefore, 



(B.15) 



where is a normahzation factor which is fixed by the operator algebra and the structure 
constants (see Table 2): 

1 



N = C4C5 



r(|)r3(|: 



Using 



with 



we find: 



2^r(|)r3(|) • 

dh l^p'^lz - Ip^ = -S{b) B{1 + a,l+b)B{l+b,-l~a- h) 
r(a)r(/5) 



S{x) = sin TTx 



22 ji 
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(B.16) 



.5/ - ^5^ 

In the second case, one screening operator is needed and the integral takes the following 
form: 



Z(a, b, c, d, e 



d'w I d'z\z?''\l - z\''\w - z\'^\w\"'\l - w 



2b I 



|2cL,,|2(i| 



|2e 



(B.17) 



The exponents a, b, c, d, e are computed in the Coulomb gas formalism |Q . Notice that 
the exponent d is in fact d' — s/2, where s comes from the Mellin regularization scheme 
(see Appendix A) . It can be shown using transformations in the complex plane that 



Z{a, b, c, d, e) 



S{a + c)~^ J d^wlwl^'^ll - wl^"" X 

S{a + b + c)S{b)\Ii{a, b, c; w)\^ + S{a)S{c)\l2{a, b, c; w)\^ 



where 



h{a,b, c; r]) 
/2(a,6, c;r]) 



2Fi(— c, —a — 6 — c — 1; —a — c; rj) 



T{-a-b-c-l)T{b+r 
r(— a — c) 

V-"^^^%^^^^^2F,(-6,a + l;a + c + 2;,) . 



r(a + c + 2) 

The generalized hypergeometric functions are defined by 

pFq{ai, ■ ■ ■ ,ap;bi, ■ ■ ,bq; z) = ^ -tttJ-. ^ 



(B.18) 



fc=0 



fc!(&l)fc ■ ■ ■ (&g)fc 



with (a)fe = r(a + A;)/r(a). 

There are several ways to compute 2' (a, 6, c, d^ e). We have used two equivalent meth- 
ods in order to have a non-trivial check. Let us explain briefly them for completeness. 
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• The first method is the one considered in |?^]. Using monodromy properties of the 
integral, we can write Z as : 



Z — Z\\j\ + Z22J2 + Z12J1J2 



with 



Ji = j z'^{l- zyh{a,h,c,z) = B{h+l,-a-h- c~l)B{d+l,e + l) 


X3-F2(— c, —a — 6 — c — —a — c, 2 + d + e; 1) , 



(B.19) 



and 



J2 = J z'^{l-zyi2{a,b,c,z) = B{a+l,c+l)B{2 + a + c + d,e + l) 


x^F2{-b, a + 1,2 + a + c + d; a + c + 2,3 + a + c + d + e;l) . (B.20) 
In these formulae the defined by: 

zii = -^S{a + cy^S'^{c + d + e)S'^{a + b + c + d + e)x 
S{b)S{a + b + c)S{d) X 

{S{b - c - d) - S{b + c - d) + S{b + c + d) - S{2a + b + c + d) 

-S{b + c + d + 2e) + S{2a + b + 3c + d + 2e)) (B.21) 

Z22 = -S{a + c)"^S'"^(c + d + e)S~'^{a + b + c + d + e)x 
S{a)S{c)S{a + c + d) x 

{S{a + b - d) + S{a - b + d) - S{a + b + d) + S{a + b + 2c + d) 

-S{a -b-d-2e) - S{a + b + 2c + d + 2e)) (B.22) 
Z12 = 2S{a + c)"^5'"^(c + d + e)S'"^(a + b + c + d + e)x 

S{a)S{b)S{c)S{a + b + c)S{d)S{a + c + d) (B.23) 



• The second method is the one considered in |57|. It consists in decomposing Z into 
its holomorphic and antiholomorphic parts: 



I = s{b)s{e) JtJi+JtJ2 + s{b)s{e + c)J:^J2 + s{b + c)s{e)J^J^ 



(B.24) 



where 



Ji = J{a, b, c, d, e) ; J2 = Jip, a, c, e, d) ; 

= J{b, —2 — a — b — c,c,e,—2 — d — e — c) ; 
J2 = J{—2 — a — b — c,b,c, —2 — d — e — c,e) , 



(B.25) 
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with the notation 

1 1 



J(a, b, c, d,e) = J duj dv u^+^+'+'il - m) V(l - vy{l - uvY = 5(1 + c, 1 + rf) x 



xB{2 + a + d + c,l + 6)3F2(-e, 2 + c + d,l + d,3 + a + b + c + d,2 + c + d;l). 



The J integrals appearing in ( |B.24| ) are not independent. Using a contour deformation it 



can be shown that they are related as 

S{a + b + c)J{ + S{a + 6) Js" = ^}''] (S{d)Jt + S{c + d)J+) ; (B.26) 

S{c + d + e)J^ + S{d + e) Jf = -wr^^ r (S{a)4 + S{a + c) J+) . (B.27) 

b[a + b + c) ^ ' 

This has the advantage that some of the Jf can be computed in an easier way than others. 

In most cases, we were able to compute exactly these types of integrals in terms of a 
product of Gamma functions by using some known relations of hypergeometric functions 



at argument z = \ [|58[ . Given their cumbersome expressions, the results are not reported 
here. When no close forms were available, the integrals were determined numerically 
(using the two different representations above) by directly calculating the hypergeometric 
functions 3-F2(-2) at z = 1. In this case, we used fast convergent expressions of the 
hypergeometric functions in order to reach rather accurate results (with approximately 
0.5 % of confidence). 

We have performed the above set of calculations for the magnetic, thermal and sub- 
magnetic perturbations. The results are in Tables |^ |TB| and |TB|. We did not pursue this 



calculation for the vacancy density perturbation for the UV difficulties explained in the 
text and also because in this case Form Factors were expected to provide a reasonable 
approximation of the correlators in all range of r. 

C Regularization of the VEV's. 

The VEV's of primary operators $/ ^ in the minimal models hAp^pi perturbed by an inte- 



grable relevant operator 0j have been conjectured in |^0|. They can be written as 

(0«|$i,fc|0,) = %,.),M(e,^)^7^ , (C.28) 

where s labels the different vacua. The prefactor -B(i,fc),j,s(^; ''?) can be further decomposed 
as 

S(i,fc),.,.(e, = bii,k),AO Q^i^, v) , (C.29) 
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where 6(i,A;),i,s(0 is a simple function for any perturbation 0j and any vacuum |0s) whereas 
QiiiiV) can be expressed as 



f+oo 
JO 



(C.30) 



In Section 4 we have presented the exphcit formulae for the three integrable deformations 
of the TIM. The above integral may diverge for some values of the two parameters 



P 



pi — p 



r^ = {^ + l)l-^k 



In fact, due to the asymptotic behavior of Fi(t; ^, rj) 



~ e 



t +00 



(C.31) 



(C.32) 



this occurs when ai{^,ri) < 0. In this case, the integral Jj(^,?7) needs to be regularized in 
order to extract its physical value. Its regularization may be performed in a number of 
equivalent ways, for instance by means of an analytic prolongation in ^ (from values for 
which Qi{^,r]) converges), possibly also using a sufficient number of times the reflection 
equations satisfied by the VEV. Here we present another simple method of controlling the 
divergences of /j(^,?7), specializing our discussion to the VEV ($2,1)4 of the sub-leading 
magnetization operator of the TIM in the massive phase reached by the perturbation of 
the vacancy operator 0i 3 = t = ip^, with < 0. This example, somehow, presents all 
possible types of divergences of the above integrals. Let us consider initially the general 
expression of the function F in the case of 0i,3 perturbation PD 



cosh2t sinh(77 — l)t 



7] 



1 



-4t 



t V2 cosht sinh(^t) sinh(^ + l)t 2^(^ +1) / 
In our example ^ = 4 and f] = Q, hence 

1 / cosh2t sinhTt 7 



(C.33) 



/4(4,6) 



dtFi{t; 4,6) 



F4(t;4,6) 



-4t 



and the asymptotic behavior 



F4(t;4,6) ~ e 



4t 



t \2cosht sinh4t 8 



t — > +00 



(C.34) 



(C.35) 



leads to the divergence of the above integral. The complete VEV under investigation can 
be expressed as 



(O.la'lO,) 



sm |7rs 
sin |s 



7 

8 


21, , 


r( 


i) 


r( 


-I) 






r 


(t 


)r 
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32 



Q4(4,6) 



(C.36) 
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where s = 1,2,3 labels the three different vacua. Assuming a regularization of Q4(4, 6), 
notice that for s = 2 the first term in ( |C.36|) implies 



(0,|a'|0,) = . (C.37) 
To compute /4(4, 6), let us break the integral as 

rl r+oo 

74(4,6) = j dtF4{t;A,6) + J dtF^{t;A,6) . (C.38) 
The first integral is always convergent since, in general 

limFi,3(t;e,^) = const , (C.39) 



thanks to a compensation between the two terms in (|C.33D . Hence it is sufficient to make 
an analytic prolongation of the second integral by subtracting its divergent part. This 
can be done by expressing initially the second integral as 

dtF,{t; 4, 6) = iF,{t; 4, 6) - — + — j + F(-4) - Y{-2) , (C.40) 

and then by making an analytic prolongation of the function 

Y{a) = / -e-'^* , 3fJe(a) > , (C.41) 
Ji t 

to the domain 3?e(a) < 0. Notice that although the analytic extension of the function 
Y(a) is not monodromic (a = is in fact a branch-cut point), its exponential e^^""^ is 
however uniquely defined and this is precisely the expression which enters the formula 
( |C36| ) for the VEV. In the punctured complex plane a G C — {0} we have 



7 



e^W = ^e-fi-) , (C.42) 

where 7 = 0.577216 ... is the Euler-Mascheroni constant and the holomorphic function 
/(a) is defined by the power series 

00 n 

/(a) = ^(-l)'^^ , aeC. (C.43) 

The last series is fastly convergent and allows good numerical estimations. 

The divergences which appear in the expression of the VEV proposed in ref. |^0|] can 
be generally tamed as in eq. (|C.40| ) above, i.e. by subtracting the leading (and sublead- 
ing) exponential terms and then performing the analytic continuation of Y{a). However, 
sometimes the first term in the r.h.s. of ( |C.40| ), i.e. the integral accompanied by the 
subtractions, is still logarithmically divergent. This is in particular the case of our exam- 
ple of {(ps)^- When this happens, the pure power law behavior in the coupling constant 
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Qi of the VEV ( |(J.28| ) gets modified. To face this situation, one may perform an extra 
regularization of the integral by making a shift of the parameter ^ which characterizes the 
minimal models 

e^e + e, (C.44) 

and then carefully taking the limit e ^ 0. In our example this results in shifting 4 — 4 + e, 
6 — > 6 + e and considering the expression 



J. 



F,it; 4 + e,6 + 6)- — + — - — dt + F (-4) - F (-2) + F (e) . (C.45) 



The limit e — > of the first integral in the previous equation is obtained by putting e = 
and hence it can be calculated by performing a numerical integration. Concerning the 
last term, once used to evaluate the VEV, we have 

e^(^) = ^ + 0{1) . (C.46) 

Let us consider now the dependence of the VEV on the coupling constant g = —g^ once 
the shift (|C.44|) has been performed. We have 

35 35 I 23_ 

g32 — >.^32"^128^ ^ (C.47) 

since the conformal weights depend on ^. By expanding this formula in powers of e we 
have 



35, J^, / N JJLf 35 / \ 35 

C, 32+ 128^ = (^O) 128'c/32 _ = ^ 32 



128 go 



(C.48) 

where we have introduced an arbitrary value of the coupling constant go for taking the log- 
arithm of the adimensional quantity ^. Once this expression is multiplied with eq. ( |C.46|) , 
it gives rise to 

In conclusion, in the example of the VEV ( |C.36| ), apart from a divergence ^ which can be 
discarded, it presents a logarithmic part and a power-law term: 

35 Q 35 

(0,|a'|0,) = 5(2,i),4,.^^ ln- + C(2,i),4,.^^. (C.50) 

9o 

It is important to notice that in the previous expression the constant -B(2,i),4,s is uniquely 
determined by this procedure and is not affected by a change of the reference coupling 
constant go- Instead, the second constant C(2,i),4,s may be freely modified by rescaling the 
arbitrary value of the reference coupling constant go- In this sense, the constant -B(2,i),4,s 
possesses a precise physical value in QFT and it may enter the definitions of amplitude 
ratios as B prefactor. 
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D Form Factors in the thermal sector 



In this appendix, we will discuss some features of the Form Factors in the thermal defor- 
mation of the TIM. In particular, we will show that it is necessary to determine indepen- 
dently (for instance by a numerical method) the one-particle Form Factors of the ipi in 
order to compute its higher particle Form Factors. We will discuss these matrix elements 
in the high-temperature phase of the model with the mass spectrum and the Z2 quantum 
number of the particles given in Table g Let 



71, 



{0\amA^,{e,...A^SOn)) . 



(D.51) 



Since (j{x) is a Z2 odd operator, the non-vanishing matrix elements will be those on Z2 
odd multi-particle states. Hence, the matrix elements which contribute to the summation 
( [4. 361 ) are, in increasing order of total energy of the corresponding states: Fi, F3, F12, -Fg, 
-^14, -Fill, F23, -F15, . . .. According to the analysis of ref. [| 
Factors can be conveniently written as 



47|, the two-particle Form 



j-rmm 



(0) 



(D.52) 



where -F^*" is the minimal solution of the set of Watson equations which has neither poles 
neither zero and and Dab are polynomials in cosh^^. The latter takes into account 
the pole structure of this matrix element (independent of the field) whereas the former 
depends on the field under consideration, in this case a{x). To determine the above 
quantities we need initially the expression of the elastic two-particle ^'-matrix of the 
model |21|, that can be expressed as 



where 



Sab{e)= n /"W"" 



/n\ — ^^^^ + ^'^'^ 



tanh 2 (^ ~ ina) 



(D.53) 



(D.54) 



The different amplitudes can be found in Table 2 of ref. and are not reported here. 
The exponents pa denote the multiplicities of the corresponding poles (located at = iira 
and 9 = in{l — a)) identified by the indices a. Correspondingly, FJ^*" is parameterized 
by: 



Tpmtn 
^ab 



(0) 



—i sinh 



n 9c.{ey 

a&Aab 



where ga{d) is given by the integral representation 



ga{9) = exp 



°° dtcoshUa - 1 2)t] , , 

■ , ,^ ■ V sin^ gt/27r 

t cosht/2sinht ^ ' ' 



(D.55) 



(D.56) 
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with 9 = m — 9. The polynomials Dab{9) are entirely determined from the poles of the 
S'-matrix. According to [|10|, they are given by 

Dab{e)= n , (D.57) 

a&Aab 

ia = n + l , ja = n , if = 2n + 1 ; 

(U.58j 

ia = n , ja = n , if pa = 2n , 
where Aab and Pa are defined in eq. ( p.53| ) and can be read from Table 2 of ref . . The 
functions 

^ cos vra — cosh 6^ 

^aie) = -- — — — , P.59 

2 cos"^ ^ 

give a suitable parameterization of the pole of the FF at 6' = in a. The asymptotic 
behavior of ga{9) and Va{9) is given by 

ga{9) ~ ; Va{9) ~ e'^' for ^ oo . (D.60) 

An upper bound on the maximal degree of the polynomial Qj^f, can be fixed by the con- 
straint []TU| 

y < A,, (D.61) 

where y is defined by 

hm F,,,...,,„(^i,...,e„)~e^l^«l . (D.62) 

\0i\^oo 

Collecting all the above results, let us consider the two-particle Form Factor ^^2(6'): 

By using eqs. (p.61|) and (|D.62| ), for the degree 6 of Q12 we have 5 < 1. The residue 
equations ( [4.42|) at the simple order poles corresponding to the bound states supply us 
with two equations, namely 

- z hm^ {9 - iuI,)F:,{9) = 7:,F; , (D.64) 

ao 

with a = 1, b = 2, c = 1, 3 and 7^^ is given by 

limJ9-zu:,)Sab{9) = {i:b? ■ (D-65) 

ao 

These two equations are able to fix unambiguously provided and are known. 
However, there is no way to determine these one-particle Form Factors in the bootstrap 
program. Notice, in fact, that the above equations are also satisfied by the two-particle FF 
of the subleading magnetization a'{x). Therefore, one needs to extract the one-particle 
FF of cr(x) and cr'(x) by means of some other independent method, as for instance the one 
provided by the numerical Truncated Conformal Space Approach, discussed in Section 5. 
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Figure 3: One-particle FF (0|(^2(0) |^i)i as a function of MR for the magnetic perturba- 
tion of the TIM. 




Figure 4: Effective LG potentials associated to the tlieory: (a) at the tricritical point 
and perturbed by (b) the leading magnetic field (gi 7^ 0), (c) the sub-leading magnetic 
field gs > 0, the leading energy density with (d) ^'2 > or (e) g2 < 0, the vacancy density 
with (f) > or (g) < 0. 
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O 0.5 1 1 -S 2 2.5 

mr 

Figure 5: The continuous line gives the UV approximation of the correlator {(p2{x)(p2{0))i 
whereas the dashed line depicts the IR approximation. An overlap of the two curves is 
observed around mr ~ 1. 
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O 0.5 1 1 .5 2 2.5 

mr 

Figure 6: The UV approximation of the correlator {(pi(x)(pi{0))2 is given by the continuous 
line whereas its IR approximation is given by the dashed line. An overlap of the curves 
is observed around mr ^ 1. 
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